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LONG-TERM  GOALS 

The  development  of  new  geoacoustic  inversion  methods,  their  use  in  the  analysis  of  shallow  water 
experimental  data,  and  evaluation  of  geoacoustic  model  and  parameter  uncertainties  including  the 
mapping  of  these  uncertainties  through  to  system  performance  uncertainties. 

OBJECTIVES 

The  development  of  methods  for  estimating  the  entire  posteriori  probability  densities  of  the 
geoacoustic  parameters  being  investigated  along  with  the  mapping  of  these  parameter  uncertainties 
through  to  characterizations  of  applied  interest  (e.g.  transmission  loss),  the  development  of  new 
geoacoustic  inversion  procedures  for  use  into  the  kHz  frequency  regime,  the  use  of  ambient  noise  for 
initial  estimation  of  seafloor  layering  structure,  and  the  demonstration  of  these  methods  in  the  analysis 
of  data  collected  during  the  Shallow  Water  2006  (SW06)  experiment. 

APPROACH 

Geoacoustic  inversion  involves  a  number  of  components:  (a)  representation  of  the  ocean  environment, 
(b)  the  inversion  procedure  selected  (e.g.  genetic  algorithm  or  simulated  annealing)  including  the 
forward  propagation  model  implemented,  and  (c)  the  estimation  of  uncertainties  associated  with  the 
parameter  estimates.  The  latter  is  critical  to  facilitate  the  mapping  of  these  uncertainties  into 
characterizations  of  applied  interest  including  the  prediction  of  total  system  performance. 

The  reporting  of  geoacoustic  parameter  estimates  without  their  associated  uncertainties  is  of  limited 
value.  Of  substantial  greater  utility  is  the  complete  a  posteriori  probability  density  (in  general,  the  joint 
density  between  all  parameters  being  estimated).  One  significant  benefit  of  obtaining  accurate  a 
posteriori  densities  of  the  geoacoustic  parameters  is  the  potential  to  map  these  through  to 
characterizations  of  applied  interest  (e.g.  transmission  loss,  source  detection  and  localization 
performance,  etc.)  in  order  to  quantify  those  uncertainties  as  well. 
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Substantial  experience  exists  in  the  application  of  full-field  geoacoustic  inversion  methods.  These  have 
been  implemented  in  a  number  of  geometries  (e.g.  fixed  vertical  and  horizontal  arrays,  towed  arrays, 
and  sonobuoys)  and  have  been  shown  to  work  well  at  low  frequencies  (<  1  kHz).  The  application  of 
these  methods  at  higher  frequencies  (into  the  few  kHz  frequency  regime)  is  at  an  early  stage.  New 
methods  are  required  which  are  robust  to  modest  geoacoustic  heterogeneity  (seafloor  parameters  as 
well  as  bathymetry)  and  temporal  fluctuations  (sound  speed  structure,  surface  waves,  and  array 
dynamics). 

Ambient  noise  provides  a  natural  illumination  source  that  can  be  used  for  waveguide  parameter 
estimation  purposes.  Our  specific  interest  here  is  in  the  use  of  ambient  noise  to  provide  initial 
estimates  of  seafloor  layering  structure.  Rough  estimates  of  seafloor  bathymetry  and  sediment 
thickness  are  needed  to  parameterize  the  waveguide  model  prior  to  carrying  out  a  geoacoustic 
inversion  procedure  regardless  of  the  type  of  data  used  for  the  inversion  itself  (e.g.  source  tow  data,  the 
radiated  signature  of  ships  of  opportunity,  or  ambient  noise). 

The  Shallow  Water  2006  experiment  took  place  in  July-September  2006  on  the  outer  edge  of  the  New 
Jersey  continental  shelf  in  approximately  80  m  deep  water.  Both  narrowband  and  broadband 
transmissions  (source  tows  and  stations)  were  made  over  a  wide  range  of  frequencies  (50  Hz  -  5  kHz) 
including  detailed  measurements  of  seafloor  structure  and  water  column  variability.  These  data  are 
available  for  geoacoustic  inversion  purposes  and  the  investigation  of  how  nuisance  parameter 
uncertainty  (e.g.  water  column  sound  speed  variability)  couples  into  seafloor  parameter  uncertainty. 

WORK  COMPLETED 

The  Shallow  Water  2006  experiment  took  place  in  July-September  2006.  One  component  of  our 
analysis  has  focused  on  short  range  broadband  geoacoustic  inversion  [1,  6,  7]  when  the  source  was 
approximately  200  m  from  the  YLA.  Due  to  the  presence  of  complex  microstructure  in  the  thermocline 
of  the  oceanic  sound-speed  profile,  fluctuations  both  in  amplitude  and  arrival  time  of  the  direct  path 
arrivals  were  observed.  Time  variation  of  the  water-column  environment  was  also  evident  during  the 
source  transmissions.  To  mitigate  the  effects  of  the  ocean  environment  on  the  seabed  property 
estimation,  a  multistage  optimization  inversion  was  employed. 

One  application  of  passive  estimation  of  the  time-domain  Green’s  function  is  in  the  use  of  cross¬ 
correlations  of  upward  and  downward  pointing  VLA  beams  observing  ambient  noise  to  extract  seabed 
layer  structure  (i.e.  a  passive  fathometer)  [2,  5].  The  use  of  adaptive  beamforming  techniques  is 
advantageous  since  at  lower  frequencies  the  horizontal  component  of  the  ambient  noise  field  can  be 
significant  [2],  This  so  called  passive  fathometer  technique  exploits  the  naturally  occurring  acoustic 
sounds  generated  on  the  sea-surface,  primarily  from  breaking  waves.  The  method  is  based  on  the 
cross-correlation  of  noise  from  the  ocean  surface  with  its  echo  from  the  seabed,  which  recovers  travel 
times  to  significant  seabed  reflectors.  To  limit  averaging  time  and  make  this  practical,  beamforming  is 
used  with  a  vertical  array  of  hydrophones  to  reduce  interference  from  horizontally  propagating  noise. 
The  initial  development  used  conventional  beamforming,  but  significant  improvements  have  been 
realized  using  adaptive  techniques.  Adaptive  methods  for  this  process  are  described  in  [2]  and  applied 
to  several  data  sets  to  demonstrate  improvements  possible  as  compared  to  conventional  processing. 

An  analytical  model  is  presented  in  [5]  for  the  passive  fathometer  response  to  ocean  surface  noise, 
interfering  discrete  noise  sources,  and  locally  uncorrelated  noise  in  an  ideal  waveguide.  The  leading 
order  term  from  the  ocean  surface  noise  produces  the  cross-correlation  of  vertical  multipaths,  yielding 
the  depth  of  sub-bottom  reflectors. 
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Finally,  we  have  explored  incorporating  Kalman  and  particle  filter  tracking  techniques  into  the 
geoacoustic  inversion  problem  [3,4],  This  enables  spatial  and  temporal  tracking  of  environmental 
parameters  and  their  underlying  probability  densities,  making  geoacoustic  tracking  a  natural  extension 
to  geoacoustic  inversion  techniques. 

RESULTS 

In  many  cases,  it  is  of  interest  to  estimate  geoacoustic  parameters  over  a  larger  spatial  region  rather 
than  just  the  parameters  characterizing  propagation  between  a  fixed  source  and  receiver  (or  receiving 
array)  location.  Data  might  be  available  at  a  moored  vertical  receiving  array  from  a  towed  acoustic 
sound  source  or  a  source  might  be  received  by  a  towed  horizontal  array.  In  both  cases,  the  typical 
approach  would  be  to  treat  each  record  of  data  independently  of  the  others  and  carry  out  a  full 
geoacoustic  inversion  for  every  record  resulting  in  a  sequence  of  geoacoustic  parameter  estimates  and, 
in  some  cases,  posteriori  probability  densities  of  the  environmental  parameters.  The  latter  enables  the 
environmental  uncertainty  to  be  projected  into  other  waveguide  characterizations  such  as  propagation 
loss  and  its  uncertainty. 

Rather  than  treating  the  data  records  independently  of  one  another,  the  geoacoustic  inversion  problem 
can  be  reformulated  as  tracking  the  evolution  of  these  parameters  and  their  associated  uncertainties  in 
space  and  time.  This  is  achieved  by  merging  geoacoustic  inversion  techniques  with  sequential  filtering 
such  as  the  Kalman  and  particle  filters,  see  Figs.  1-2.  The  interaction  between  the  environmental 
parameters  and  the  acoustic  field  can  involve  a  high  level  of  nonlinearity.  In  addition,  the  posteriori 
probability  densities  of  geoacoustic  parameters  can  be  non-Gaussian.  Thus,  geoacoustic  tracking 
requires  tracking  filters  that  can  handle  nonlinear,  non-Gaussian  systems.  We  have  demonstrated  the 
suitability  of  three  such  filters  in  geoacoustic  tracking-  the  extended  Kalman  filter  (EKF),  the 
unscented  Kalman  filter  (UKF),  and  the  particle  filter  (PF)  on  real  data  from  the  SWellEx-96 
experiment  [3]. 

In  a  review  paper  we  have  studied  the  basis  and  use  of  sequential  filtering  in  ocean  acoustics  [4], 
Sequential  filtering  provides  an  optimal  framework  for  estimating  and  updating  the  unknown 
parameters  of  a  system  as  data  become  available.  Despite  significant  progress  in  the  general  theory  and 
implementation,  sequential  Bayesian  filters  have  been  sparsely  applied  to  ocean  acoustics.  The 
foundations  of  sequential  Bayesian  filtering  with  emphasis  on  practical  issues  are  first  presented 
covering  both  Kalman  and  particle  filter  approaches.  Filtering  becomes  a  powerful  estimation  tool, 
employing  prediction  from  previous  estimates  and  updates  stemming  from  physical  and  statistical 
models  that  relate  acoustic  measurements  to  the  unknown  parameters.  Ocean  acoustic  applications  are 
then  reviewed  focusing  on  the  estimation  of  environmental  parameters  evolving  in  time  or  space.  Some 
possible  scenarios  for  geoacoustic  inversion  are  shown  in  Fig.  3. 

We  demonstrated  the  use  of  PF  on  event  S9  from  SWellEx-96  [3],  The  environment  is  given  in  Fig.  4. 
Event  S9  is  selected  here  since  the  track  is  perpendicular  to  bathymetric  lines  giving  the  highest  rate  of 
change  of  the  environment  as  the  source  ship  moved  into  shallower  water.  Even  though  the  source  is 
successfully  tracked  by  using  a  200-particle  PF,  a  second  filtering  is  performed  using  a  10  000-particle 
PF.  This  enables  us  to  obtain  much  better  histograms  of  particles  that  represent  the  marginal  posterior 
PDFs  of  the  geoacoustic  parameters.  Evolving  marginal  posterior  PDFs  of  the  four  geoacoustic 
parameters  are  given  in  Fig.  5.  The  water  depth  is  well  determined  throughout  the  track  and  the 
posterior  PDFs  correspond  to  narrow  Gaussian  PDFs.  The  slope  of  the  sediment  sound  speed  profile  is 
likewise  well-determined,  closely  following  the  underlying  true  values  obtained  from  a  SWellEx 
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region  database.  Note,  however,  that  both  the  sediment  thickness  and  top  layer  sound  speed  have 
highly  non-Gaussian,  quickly-evolving  posterior  PDFs  with  multiple  peaks. 

IMPACT  /  APPLICATIONS 

Geoacoustic  inversion  techniques  are  of  general  interest  for  the  estimation  of  waveguide  parameters 
thus  facilitating  system  performance  prediction  in  shallow  water.  Natural  transition  paths  for  these 
results  will  be  the  PEO-C4I  Battlespace  Awareness  and  Information  Operations  Program  Office 
(PMW-120)  and  the  Naval  Oceanographic  Office. 

RELATED  PROJECTS 

This  project  is  one  of  several  sponsored  by  ONR  Code  3210A  to  participate  in  the  Shallow  Water 
2006  experiment  and  participate  in  the  analysis  of  the  resulting  data. 
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Figure  1.  Sequential  Bayesian  filtering.  From  state  Xk-i,  state  x*  is  first  predicted  via  the  state 
equation,  providing  x*| t-i  .  As  data  yk  becomes  available,  the  observation  equation  is  employed  to 

update  state  Xk\k-i ,  providing  Xk\k  • 


Figure  2.  A  quick  guide  to  filter  selection  leading  to  the  Kalman  filter  (KF),  extended  Kalman  filter 
(EKF),  unscented  Kalman  filter  ( UKF '),  and  particle  filter  (PF). 
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Figure  3.  Geoacoustic  environmental  tracking:  (a)  Temporal  tracking  of  the  ocean  sound  speed 
profile  for  a  fixed-receiver  and  a  fixed-source  and  (b)  tracking  of  the  changing  environment 
between  the  receiver  and  a  moving  source.  Here  shown  for  a  vertical  line  array  (VLA)  of  receivers. 


Figure  4.  The  environmental  parameters  used  to  model  the  vertical  profiles  at  the  VLA  and  source 
for  event  S9  from  SWellEx-96.  The  values  in  parentheses  are  at  the  location  of  the  VLA. 
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Figure  5.  Evolving  marginal  posterior  PDFs  for  the  geoacoustic  parameters:  (a)  water  depth  (wd), 
(b)  sediment  thickness  (hsed),  (c)  sediment  sound  speed  at  the  water-sediment  interface  (cs),  and  (d) 
sediment  sound  speed  gradient  (dcfck,).  Vertical  slices  of  the  PDFs  are  given  at  5  min  intervals. 

Solid  lines  represent  the  true  values. 
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